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METHOD, APPARATUS AND PROGRAM 
FOR RESTORING PHASE INFORMATION 

5 

BACKGROUND OF THE INVENTION 
Field of the invention 

The present invention relates to a method, an apparatus , 
and a program for restoring phase information to be used for 

10 constructing an image on the basis of image information 
obtained by radiation imaging . In this application , the word 
"radiation" is used in a broad sense that includes corpuscular 
beams such as an electron beam, and electromagnetic waves, 
in addition to X-rays, a-rays, [i-rays, y-rays, ultraviolet 

15 rays and the like. 



Description of a Related Art 

Conventionally, an imaging method using an X-ray, etc. 
is utilized in various fields , and particularly in the medical 

20 field, the method is one of the most important means for 
diagnosis . Since the first X-ray photograph was implemented, 
the X-ray photography method has been repeatedly improved, 
and a method using a combination of a fluorescent screen and 
an X-ray film is in the mainstream at present. On the other 

25 hand, in recent years, various digitized devices applying 
X-ray CT, ultrasonic waves, MRI , etc. are in practical use, 
and construction of a diagnostic information processing 
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system or the like in hospitals is being promoted. With 

respect to X-ray images, many studies are also made for 

digitizing the imaging system . Digitizing the imaging system 

enables to store a large amount of data for a long period 

5 without degradation in image quality, and helps to make 

progress toward the medical diagnostic information processing 
system. 

Now, a radiation image obtained as described above is 
generated by converting intensity of the radiation 
10 transmitted through an object into brightness of the image. 
For example, when performing imaging on a region including 
a bone part, the radiation transmitted through the bone part 
is largely attenuated, and the radiation transmitted through 
a region other than the bone part, i . e . , a soft part is slightly 
15 attenuated. In this case, since the difference in intensity 
between the radiations transmitted through different tissues 
is large, a high-contrast radiation image can be obtained. 

On the other hand, for example, when imaging a region 
of the soft part such as a breast, since radiation is easily 
20 transmitted through the soft part as a whole, the difference 
between tissues in the soft part is difficult to appear as 
the difference in intensity of the transmitted radiation. 
Therefore , only a low-contrast radiation image can be obtained 
with respect to the soft part. Thus, the conventional 
25 radiation imaging method is not suitable as a method of 
visualizing the slight difference between tissues in the soft 
part. 
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Here, as information included in the radiation 
transmitted through the object, there is phase information 
in addition to intensity information. Recently, a 
phase-contrast method of generating an image using this phase 
information is under study. The phase-contrast method is an 
image construction technology to convert phase difference 
produced by an X-ray, etc. transmitted through the object 
into brightness of an image . 

The phase-contrast method includes a technique for 
obtaining phase difference on the basis of interference X-ray 
produced by using an interferometer or a zone plate, and a 
technique for obtaining phase difference on the basis of 
diffraction X-ray . Among these techniques , the technique for 
obtaining phase difference on the basis of diffraction X-ray 
(diffraction technique) is to obtain phase difference on the 
basis of the following principle. An X-ray, for example, 
propagates within a material because a wave progresses 
similarly to light. The propagation rate varies according 
to the refraction index of the material. Therefore, when 
applying an X-ray having uniform phase toward the object, 
the propagation rate of the X-ray varies according to the 
difference between tissues in the object. Thereby, the 
wavefront of the X-ray transmitted through the object is 
distorted and, as a result, diffraction fringes are produced 
in the X-ray image obtained on the basis of the transmitted 
X-ray. The pattern of the diffraction fringes differs in 
accordance with the distance between the screen on which the 




X-ray image is formed and the object, and the wavelength of 
the X-ray. Thus, analyzing two or more X-ray images having 
different patterns of diffraction fringes , phase differences 
of the X-ray generated in the respective positions on the 
5 screen can be obtained. By converting the phase differences 
into brightness, an X-ray image that clearly shows the 
difference between tissues in the object can be obtained. 

Particularly, in the radiation after transmitted through 
a soft part of the object, since the difference in phase is 
10 larger than the difference in intensity in accordance with 
the difference between tissues through which the radiation 
has been transmitted, the subtle difference between tissues 
canbe visualizedby using the phase-contrast method. In order 
to use the above-described phase-contrast method, imaging 
15 conditions in radiation imaging or techniques for restoring 
phase information from patterns of diffraction fringes are 
under study. 

B. E. Allman et al . "Noninterferometric quantitative 
phase imaging with soft x rays" , J. Optical Society of America 

20 A, Vol. 17, No. 10 (October 2000), pp. 1732-1743 discloses 
that an X-ray image is constructed by performing phase 
restoration on the basis of image information obtained by 
imaging with soft X-rays. In this reference, the basic 
equation of phase restoration, TIE (transport of intensity 

25 equation) is used. 

dl(r) 

* dz --V x -{I(r) V x <ft (r)} •••(!) 
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Where r = (x, y, z) is vector, and 

In addition, k denotes wave number. 

Next, a principle of the phase restoration will be 
described by referring to Fig. 21. As shown in Fig. 21, an 
X-ray having a wavelength X is input from the left side of 
the drawing, transmitted through an object plane 1 and enters 
a screen 2 at a distance of z from the object plane 1. Here, 
assuming that the intensity and the phase of the X-ray in 
a position (x, y) on the screen 2 are I(x, y) and <|>(x, y) , 
respectively. In this case, relationship expressed by the 
following equation holds between the intensity I(x, y) and 
phase <|)(x, y) . Here, the intensity I is square of amplitude 
of wave. 

27t dl<x,y) 

15 "A Jl V ■<!(*,*> V<Mx,y)} .-.(2) 

Substituting k = 2tiA into Eq. (2) and rewriting (x, y) 
components into vector r, TIE expressed by Eq. (1) is derived. 

However, since the above TIE is difficult to be solved, 
TIE has been used mainly by performing approximation thereon . 

20 For example, T. E. Gureyev et al. "Hard X-ray quantitative 
non-interferometric phase-contrast imaging" , SPIE, Vol. 3659 
(1999), pp. 356-364 discloses that an X-ray image is 
constructed by performing phase restoration on the basis of 
image information obtained by imaging with hard X-rays . In 

25 this reference, TIE expressed by Eq. (1) is approximated as 



5 



10 



follows. First, Eq. (1) is developed. In the following 
equations, the vector r in the above reference is rewritten 
into (x, y) components. 

dx l I(x ' y) ^-J + -a7l I(x ' y) — ay— j 

= I(x,y) [ 92<t> (X,Y) | d2<t> (X,Y) ) | dKx^jrt a <j> (x,y) 
V dx 2 dy 2 J dx dx 

dl(x,Y) d 4> (x,y) 
dy dy 

= I(x,y) V 2 <t> (x,y)+Vl(x,y) • V <*> (x,y) •••(3) 
Where 



V 2 = ^— + 



5x 2 dy 2 

Approximating the second term on the right side of Eq. 
(3) to zero, the approximation equation expressed by the 
following equation (4) is obtained. 

dl(x,Y) ... I(x,y) 

Jl = " V <Mx,y) ■••(4) 



In Eq. (4) , <|>(x, y) can be obtained from I (x, y) by a solution 
method such as the finite element method. 

In addition, T. E. Gureyev et al . "Quantitative In-Line 
Phase-Contrast Imaging with Multienergy X Rays", Physical 
15 Review Letter, Vol. 86, No. 25 (2001) , PP. 5827-5830 discloses 
that an X-ray imaging is performed by using three kinds of 
X-rays having different wavelengths and the phase restoration 
is performed on the basis of the obtained image information . 
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In this reference, attention is given to the relationship 
between phase and intensity of the X-ray just after transmitted 
through the object and intensity of the X-ray in a position 
at a predetermined distance from the object . When performing 
the X-ray imaging, the structure shown in Fig. 22 is assumed. 
That is, as shown in Fig. 22, the three kinds of X-rays having 
wavelengths A 0 , X x , and A 2 , respectively, are transmitted 
through the object 3, and enter the screen 2 that is disposed 
in a position at a distance R from the object plane 1. 

In this case, assuming that r x = (x, y) , the following 
relationship is held between the intensity I(r x , 0, A 0 ) and 
the phase <J»(r ± , 0, A 0 ) of the X-ray of wavelength Ao, which 
is just after transmitted through the object 3, and the 
intensity I (r_i_, R, A^) of the diffraction X-ray of wavelength 
Am which is detected on the screen 2 . In the following Equation 
(5), I(r x , 0, A 0 ) = exp {-M(r x , 0, A 0 ) } . 

/M(r ± ,0,A 0 ) \ 
A I -V 2 * (r x ,0, A 0 ) 

\ VM-V<Mr ± ,0, A 0 ) / 

Where, 
A = | -a 




(5) 



o _ = 




, g m = ln{I(r ± ,R, A m ) } (m = 0,1,2) 



In Eq. (5) , if VX4 - V<J)(r x , R, A 0 ) is sufficiently small, 
approximation can be as follows. 
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-1 T 0 W M(r ± ,0,A 0 ) \ /g Q 

.-"x 3 °ir x J [ -V 2 *(r ± ,0,A 0 ) )-[ gi 1 •••< 6 > 

Further, from Eq. (6) , the intensity and the phase of 
the X-ray just after transmitted through the object 3 are 
expressed as follows. 

M(r ± ,0, A 0 )=-^-(g 0 -<7- 2 gi ) ... (7 ) 

-V 2 *(r ± ,0, A 0 ) = — _(ag 0 -a- 2 gi ) ... (8 ) 

Where, AA = Ai-A 0 and a = a x = Ai/Ao. 

The phase <t>(r x , R, A 0 ) can be obtained by performing 
inverse Laplacian computation on the Laplacian V 2 <j>(rj., R, A 0 ) 
of the phase in Eq . ( 8 ) . Further , a visible image representing 
the object can be obtained by converting this phase into 
brightness in the image. Thus, by using Eq. (8), the 
computation for phase restoration can be easily performed 
on the basis of a small number of irradiation images obtained 
by changing the wavelengths. Therefore, in T. E. Gureyev et 
al. "Quantitative In-Line Phase-Contrast Imaging with 
Multienergy X Rays", Physical Review Letter, Vol. 86, No. 
25 (2001) , pp. 5827-5830, X-ray imaging is performed by using 
X-rays having three kinds of wavelengths (energy) A 0 = 3.8 
A (E 0 = 3.3 keV) , A x = 7.3 A (E x = 1.7 keV) , and A 2 = 2.5 A 
(E 2 = 5.0 keV) . 

Although a thin ob j ect can be imaged at these wavelengths , 
there is a problem that the imaging method is unsuitable in 
the case where an object having a larger thickness such as 
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a breast or a chest of a human body is to be imaged because 
X-ray absorption when transmitted through the object is too 
large at these wavelengths. 

Further, in order to perform higher precision phase 
restoration, it is desired to use a high-definition (high 
resolution) screen constituted by a many number of detecting 
elements the size of which is made as small as possible. 

However, in the case where the high-definition screen 
is used, there is a problem that the X-ray imaging becomes 
vulnerable to the influence of noise . In order to reduce the 
influence of noise, it is conceivable that the irradiation 
amount of X-ray is increased . However , in the case of a living 
organism such as a human body, there is a problem that the 
irradiation amount can not be increased for fear of exposure . 
By the way, in the above -de scribed papers of B. E. Allman 
et al . and T . E . Gureyev et al . in which a non-living organism 
is imaged, the relationship between noise and irradiation 
amount is not described. 



SUMMARY OF THE INVENTION 

The present invention has been achieved in view of the 
above-described problems. A first object of the present 
invention is to provide a phase information restoring method 
that enables to improve estimation accuracy of phase by using 
radiation with energy of high transmittance when constructing 
a radiation image of a living organism such as a human body 
by the phase-contrast method. Additionally, a second object 



of the present invention is to provide a phase information 
restoring method in which the influence of noise can be reduced 
without increasing irradiation amount of X-rays . Further a 
third object of the present invention is to provide a phase 
information restoring apparatus and a phase information 
restoring program for using these phase information restoring 
methods . 

In order to solve the above described problems , a phase 
information restoring method according to a first aspect of 
the present invention is a method of restoring phase 
information on radiations transmitted through an object on 
the basis of detection data obtained by detecting intensity 
of the radiations transmitted through the object . The method 
comprises the steps of: (a) obtaining plural sets of detection 
data respectively representing plural kinds of radiation 
image information on a detection plane at a predetermined 
distance from the object by using plural radiations having 
different wavelengths with energy from 16 keV to 30 keV to 
detect intensity of the plural radiations on the detection 
plane; and (b) restoring phase information on the radiations 
transmitted through the object on the basis of the plural 
sets of detection data so as to obtain phase data. 

Further, a phase information restoring method according 
to a second aspect of the present invention is a method of 
restoring phase information on a radiation transmitted 
through an object on the basis of detection data obtained 
by detecting intensity of the radiation transmitted through 
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the object. The method comprises the steps of: (a) obtaining 
plural sets of detection data respectively representing 
plural kinds of radiation image information on plural 
detection planes at different distances from the object by 
using a radiation having a predetermined wavelength with 
energy from 16 keV to 30 keV to detect intensity of the radiation 
on the plural detection planes; and (b) restoring phase 
information on the radiation transmitted through the object 
on the basis of the plural sets of detection data so as to 
obtain phase data. 

Furthermore, a phase information restoring method 
according to a third aspect of the present invention is a 
method of restoring phase information on a radiation 
transmitted through an object on the basis of detection data 
obtained by detecting intensity of the radiation emitted from 
a radiation source having a focal spot size a in a position 
at a distance R from the object and transmitted through the 
object. The method comprises the steps of: (a) obtaining 
plural sets of detection data respectively representing 
plural kinds of radiation image information on plural 
detection planes at different distances z ± from the object 
by using a radiation detector for detecting intensity of 
applied radiation on the plural detection planes to generate 
a detection signal representing radiation image information 
in which a pixel size is not less than 7iaz/3R, where z is 
a maximum value of z i; and (b) restoring phase information 
on the radiation transmitted through the object on the basis 
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of the plural sets of detection data so as to obtain phase 
data. 

Moreover, aphase information res toringmethod according 
to a fourth aspect of the present invention is a method of 
restoring phase information on a radiation transmitted 
through an object on the basis of detection data obtained 
by detecting intensity of the radiation emitted from a 
radiation source having a focal spot size a in a position 
at a distance R from the object and transmitted through the 
object. The method comprises the steps of: (a) acquiring 
plural sets of first detection data respectively representing 
plural kinds of radiation image information on plural 
detection planes at different distances z ± from the object 
and obtained by detecting intensity of the radiation on the 
plural detection planes ; (b) respectively generating plural 
sets of second detection data by suppressing spatial frequency 
components larger than 3R/2tcctz with respect to the plural 
sets of first detection data, where z is a maximum value of 
z ± ; and (c) restoring phase information on the radiation 
transmitted through the object on the basis of the plural 
sets of second detection data so as to obtain phase data. 

A phase information restoring apparatus according to 
a first aspect of the present invention is an apparatus for 
restoringphase information on radiations transmitted through 
an object on the basis of detection data obtained by detecting 
intensity of the radiations transmitted through the object. 
The apparatus comprises : a radiation source for emitting each 
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of plural radiations having different wavelengths with energy 
from 16 keV to 30 keV; detecting means for detecting intensity 
of the radiation emitted from the radiation source and 
transmitted through the object so as to obtain detection data 
5 representing radiation image information; and phase data 
calculating means for restoring phase information on the 
radiations having different wavelengths and transmitted 
through the object on the basis of plural sets of detection 
data obtained by detecting intensity of the plural radiations 
10 so as to obtain phase data. 

Further, a phase information restoring apparatus 
according to a second aspect of the present invention is an 
apparatus for restoring phase information on a radiation 
transmitted through an object on the basis of detection data 
15 obtained by detecting intensity of the radiation transmitted 
through the object. The apparatus comprises: a radiation 
source for emitting a radiation having a predetermined 
wavelength with energy from 16 keV to 30 keV; detecting means 
for detecting intensity of the radiation emitted from the 
20 radiation source and transmitted through the object so as 
to obtain detection data representing radiation image 
information; driving means to be used for changing a distance 
between the object and the detecting means; and phase data 
calculating means for restoring phase information on the 
25 radiation transmitted through the object on thebasis of plural 
sets of detection data obtained by detecting intensity of 
the radiation at different distances so as to obtain phase 
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data. 

Furthermore, a phase information restoring apparatus 
according to a third aspect of the present invention is an 
apparatus for restoring phase information on a radiation 
transmitted through an object on the basis of detection data 
obtained by detecting intensity of the radiation, which has 
been emitted from a radiation source having a focal spot size 
a in a position at a distance R from the object and transmitted 
through the object, on plural detection planes at different 
distances z ± from the object. The apparatus comprises: 
detecting means for detecting intensity of applied radiation 
so as to obtain detection data representing radiation image 
information in which a pixel size is not less than 7iaz/3R, 
where z is a maximum value of z i; and phase data calculating 
means for restoring phase information on the radiation 
transmitted through the object on the basis of plural sets 
of detection data obtained by detecting intensity of the 
radiation transmitted at different distances z ± so as to obtain 
phase data. 

Moreover, a phase information restoring apparatus 
according to a fourth aspect of the present invention is an 
apparatus for restoring phase information on a radiation 
transmitted through an object on the basis of detection data 
obtained by detecting intensity of the radiation, which has 
been emitted from a radiation source having a focal spot size 
a in a position at a distance R from the object and transmitted 
through the object, on plural detection planes at different 
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distances z ± from the object . The apparatus comprises : signal 
processing means for respectively generating plural sets of 
second detection data by suppressing spatial frequency 
components larger than 3R/2tkjz with respect to plural sets 
of first detection data obtained by detecting intensity of 
the radiation at different distances Zi , where 2 is a maximum 
value of Zi ; and phase data calculating means for restoring 
phase information on the radiation transmitted through the 
object on the basis of the plural sets of second detection 
data generated by the signal processing means so as to obtain 
phase data. 

A phase information restoring program according to a 
first aspect of the present invention is aprogram for restoring 
phase information on radiations transmitted through an object 
on the basis of detection data obtained by emitting the 
radiations from a radiation source and detecting intensity 
of the radiation transmitted through the object . The program 
actuates a CPU to execute the procedures of: controlling the 
radiation source to emit each of radiations having different 
wavelengths with energy from 16 keV to 30 keV; obtaining a 
Laplacian of phase on the basis of plural sets of detection 
data obtained by detecting intensity of the radiations having 
different wavelengths; and obtaining phase data of the 
radiation by performing inverse Laplacian computation on the 
Laplacian of phase. 

Further , aphase information restoring program according 
to a second aspect of the present invention is a program for 
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restoring phase information on a radiation transmitted 
through an object on the basis of detection data obtained 
by emitting the radiation from a radiation source and detecting 
intensity of the radiation transmitted through the object. 
The program actuates a CPU to execute the procedures of: 
controlling the radiation source to emit a radiation having 
a predetermined wavelength with energy from 16 keV to 30 keV; 
obtaining a Laplacian of phase on the basis of plural sets 
of detection data obtained by detecting intensity of the 
radiation at different distances; and obtaining phase data 
of the radiation by performing inverse Laplacian computation 
on the Laplacian of phase. 

Furthermore, a phase information restoring program 
according to a third aspect of the present invention is a 
program for restoring phase information on a radiation 
transmitted through an object on the basis of detection data 
obtained by detecting intensity of the radiation emitted from 
a radiation source having a focal spot size ct in a position 
at a distance R from the object and transmitted through the 
ob j ect . The program actuates a CPU to execute the procedures 
of: obtaining plural sets of detection data respectively 
representing plural kinds of radiation image information on 
plural detection planes at different distances z ± from the 
object by using a radiation detector for detecting intensity 
of applied radiation on the plural detection planes to generate 
a detection signal representing radiation image information 
in which a pixel size is not less than 7iaz/3R, where z is 
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a maximum value of z ± ; and obtaining a Laplacian of phase on 
the basis of the plural sets of detection data; and obtaining 
phase data of the radiation by performing inverse Laplacian 
computation on the Laplacian of phase. 

Moreover, a phase information restoring program 
according to a fourth aspect of the present invention is a 
program for restoring phase information on a radiation 
transmitted through an object on the basis of detection data 
obtained by detecting intensity of the radiation emitted from 
a radiation source having a focal spot size a in a position 
at a distance R from the object and transmitted through the 
object . The program actuates a CPU to execute the procedures 
of : acquiring plural sets of first detection data obtained 
by detecting intensity of the radiation on plural detection 
planes at different distances Zi from the object; respectively 
generating plural sets of second detection data by suppressing 
spatial frequency components larger than 3R/2tcctz with respect 
to the plural sets of first detection data , where z is a maximum 
value of z ± ; obtaining a Laplacian of phase on the basis of 
the plural sets of second detection data; and obtaining phase 
data of the radiation by performing inverse Laplacian 
computation on the Laplacian of phase. 

According to the first and second aspects of the present 
invention, when constructing the radiation image of a living 
organism such as a human body by the phase-contrast method, 
the estimation accuracy of phase can be improved by using 
radiation with energy of high transmittance . Further, 
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according to the third and fourth aspects of the present 
invention, the influence of noise can be reduced without 
increasing irradiation amount of the X-ray by suppressing 
the predetermined spatial frequency component. 

BRIEF DESCRIPTION OF THE DRAWINGS 

Fig. 1 is a block diagram showing a constitution of a 
phase information restoring apparatus according to a first 
embodiment of the present invention; 

Fig . 2 is a diagram showing a constitution of an imaging 
unit shown in Fig. 1; 

Fig. 3 is a graphical representation showing 
transmittance as a function of X-ray energy; 

Fig. 4 is a graphical representation showing 
relationship between 8d and energy ; 

Fig. 5 is a flowchart showing a phase information 
restoring method according to the first embodiment of the 
present invention; 

Fig. 6 is a diagram showing the imaging unit being 
generating each of two kinds of X-rays having different 
wavelengths X Q and X x at the same imaging distance; 

Fi g- 7 is a block diagram showing a modified example 
of the phase information restoring apparatus according to 
the first embodiment of the present invention; 

Fig. 8 is a diagram showing a constitution of a reading 
unit shown in Fig. 7; 

Fig. 9 is a block diagram showing a constitution of a 
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phase information restoring apparatus according to a second 
embodiment of the present invention; 

Fig . 10 is a diagram showing a constitution of the imaging 
unit shown in Fig. 9; 

Fig. 11 is a flowchart showing a phase information 
restoring method according to the second embodiment of the 
present invention; 

Fig. 12 is a block diagram showing a constitution of 
a phase information restoring apparatus according to a third 
embodiment of the present invention; 

Fig . 13 is a diagram showing a constitution of an imaging 
unit shown in Fig. 12; 

Fig . 14 is a diagram for explanation of a detection signal 
obtained by detecting an X-ray generated from an X-ray source 
having a finite focus size by using detecting elements included 
in a sensor; 

Fig. 15 is a graphical representation showing a blur 
function as a function of spatial frequency; 

Fig. 16 is a flowchart showing a phase information 
restoring method according to the third embodiment of the 
present invention; 

Fi 9 - 17 is a diagram showing a constitution of the imaging 
unit when performing enlarging radiography; 

Fig. 18 is a block diagram showing a constitution of 
a phase information restoring apparatus according to a fourth 
embodiment of the present invention; 

Fig. 19 is a flowchart showing a phase information 
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restoring method according to the fourth embodiment of the 
present invention; 

Fig. 20 is a block diagram showing a constitution of 
a modified example of the phase information restoring 
apparatus according to the fourth embodiment of the present 
invention ; 

Fig. 21 is a diagram for explanation of a principle of 
phase restoration; and 

Fig. 22 is a diagram for explanation of a principle of 
phase restoration. 

DESCRIPTION OF THE PREFERRED EMBODIMENTS 

Now, referring to the drawings, embodiments of the 
present invention will be described in detail. The same 
component elements will be given with the same reference 
numerals and the descriptions thereof will be omitted. 

Fig. 1 is a block diagram showing a constitution of a 
phase information restoring apparatus according to a first 
embodiment of the present invention . As shown in Fig . 1 , this 
phase information restoring apparatus has an imaging unit 
100 for outputting detection data that represents radiation 
image information on an objectby applying X-rays to the object, 
an image constructing unit 200 for generating image data by 
restoring phase information on the basis of the detection 
data, a display unit 300 for displaying a visible image on 
the basis of the image data, and an output unit 400 for printing 
out the visible image on a film, etc. 
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Fi S- 2 is a diagram showing a constitution of the imaging 
unit 100 . The imaging unit 100 includes an X-ray source 101 , 
a monochromator 102, and a sensor 103. 

As the X-ray source, it is desirable to use one capable 
of generating a beam that is highly coherent andmonochromatic . 
Here, the highly monochromatic beam indicates a beam that 
mainly has single wavelength . However , it is unnecessary that 
the beam strictly has single wavelength. For this purpose, 
in the embodiment , the X-ray source 101 including a synchrotron 
radiation source for generating synchrotron radiation and 
the monochromator 1 02 i s used . By extracting a predetermined 
wavelength component of an X-ray included in the synchrotron 
radiation by using the monochromator 102, a monochromatic 
X-ray is obtained. In the embodiment, as described later, 
two kinds of X-rays having wavelengths X 0 and A* respectively 
are obtained. The synchrotron radiation represents 
radiation (electromagnetic wave) that is generated by 
circularly or spirally moving an electron in magnetic field, 
and has wide range continuous spectrum including an X-ray 
region . In the synchrotron radiation source generating such 
synchrotron radiation, a wavelength component included in 
the synchrotron radiation can be changed by varying 
centripetal acceleration of an electron. 

As shown in Fig. 2, the X-ray generated from the X-ray 
source 101 is transmitted through the object 10 and enter 
the sensor 103 to produce diffraction fringes . Hereinafter , 
a distance between the object 10 and the sensor 103 is referred 
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to as "imaging distance" . 

The sensor 103 is used as a screen for causing the X-ray 
to enter and produce diffraction fringes, and outputs 
detection signals representing intensity of diffraction X-ray 
entering the respective positions of the sensor 103 . As the 
sensor 103, a two-dimensional sensor such as a CCD (charge 
coupled device) is used, which has a plurality of detecting 
elements that convert intensity of the incident X-rays into 
electric signals and output the signals. 

Further, the imaging unit 100 includes an amplifier 104 
and an A/D converter 105. The amplifier 104 amplifies the 
detection signal outputted from the sensor 103. The A/D 
converter 105 converts the detection signal amplified by the 
amplifier 104 into a digital signal (referred to as "image 
signal" or "detection data") and outputs it to the image 
constructing unit 200. 

Referring to Fig. 1 again, the image constructing unit 
200 has a storage unit 201 for temporarily storing the detection 
data outputted from the imaging unit 100, a Laplacian 
processing unit 202 for calculating a value that corresponds 
to a Laplacian of phase on the basis of the two pieces of 
detection data obtained by the two kinds of X-rays having 
different wavelengths at the same imaging distance , an inverse 
Laplacian processing unit 203 for performing inverse 
Laplacian computation for phase restoration, an image 
processing unit 204 for generating image data on the basis 
of the restored phase information, and a control unit 205 
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for controlling the respective units 201-204 and wavelengths 
of the X-rays in the imaging unit 100 . The image constructing 
unit 200 may be configured with a digital circuit or software 
and a CPU. In the latter case, the control unit 205 including 
the CPU processes the detection data on the basis of a phase 
information restoring program recorded on a recording medium 
206. As the recording medium 206, a flexible disk, a hard 
disk, an MO, an MT, a RAM, a CD-ROM, a DVD-ROM, etc. are 
applicable . 

The display unit 300 is a display device such as a CRT, 
for example, and displays a visible image on the basis of 
the image data that represents the phase information restored 
by the image constructing unit 200 . Further, the output unit 
400 is a laser printer , for example, and prints out the visible 
image on a film, etc. on the basis of the image data. 

Next, the relationship between the wavelength (energy) 

of the X-ray generated in the X-ray source 101 and transmi ttance 

in the object 10 will be described. Here, the object 10 is 

assumed to have a thickness d, as shown in Fig. 2 . Expressing 

the thickness of the object as d, and refraction index as 

n = 1-8-10 (i is an imaginary unit) , the X-ray <|> OUT just after 

transmitted through the object is expressed as follows by 

using the X-ray fa just before transmitted through the object 
10. 

d> =6 e " iknd - ch ~k£d -ik(l-(5 )d /AV 

V OUT ^ IN e —(P IN e e **"(9) 

Where, assuming the wavelength of the X-ray as X, an equation 
k = 2n/X stands up. 
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Therefore, the intensity I OUT of the X-ray just after 
transmitted through the object 10 is expressed as follows 
by using the of the X-ray just before transmitted through 
the object 10. 

I out - I *out I 2 = Uxn I V 2k * d - I IN e~ 2k * d (10) 
Accordingly, the transmittance T is expressed as follows. 

T = Iqut _ -2k£d 



= e ... (ii) 



Next, as the object 10, abreast is assumed. Fig. 3 shows 
transmittance as a function of X-ray energy. Note that P is 
obtained by assuming standard fat composition and the 
thickness of the object 10 is assumed as d = 5 cm. As shown 
in Fig. 3, the lower the energy, the more drastically the 
transmittance is reduced. On this account, the X-ray energy 
(E s 3.3 keV, 1.7 KeV, and 5.0 keV) as described in T. E. 
Gureyev et al . "Quantitative In-Line Phase-Contrast Imaging 
with Multienergy X Rays", Physical Review Letter, Vol. 86, 
No. 25 (2001), pp. 5827-5830, is not suitable in the case 
where the object having a larger thickness such as a breast 
is subjected to imaging. Further, in order to improve energy 
efficiency with high transmittance, it is desired that the 
X-ray energy be on the order of 22 keV at which the transmittance 
is about 10 %. However, in order to allow the X-ray to reach 
the sensor 103 while suppressing the amount of X-ray absorption 
(exposure amount) by the object 10, equal or more than 16 
keV of X-ray energy at which the transmittance is at least 
equal or more than 1 % is required. 
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On the other hand , the phase-contrast method is an imaging 
method of detecting a phase change by separating the object 
10 and the sensor 103 in Fig. 2 at a predetermined distance. 
The phase change of the X-ray by the object 10 is expressed 
as e" lk(1 " 8)d in Eq. (8) . Here, with reference to the e" ikd of 
the phase change of the X-ray in the case where the object 
is vacuum (8=0), the phase change of the X-ray by the object 
10 is expressed by e ikSd . The phase is kSd on the shoulder 
of the exponential . 

Here, 8 has dependence on X-ray energy. Fig. 4 shows 
the relationship between 8d and energy in the case where the 
standard fat composition is assumed. Note that d = 5 cm. 

By the way, considering the TIE as a basis when phase 
estimation shown in Fig. 4 in one-dimension for simplicity, 
neglecting signs, and further arranging by substituting the 
phase $ = k8d, the TIE is expressed as the follows. 

1 dl ^ d 2 dd 
' T TT = "TT ••(12) 

I d Z d 

Here, the left side is a value representing how much 
the intensity I of the X-ray changes according to the distance 
z as a ratio relative to the intensity I , and the value obtained 
by experiments is atmost 1% (|d/l) x (dl/dz) \>10 2 ) . Further, 
the right side is a value dependent on how much the value 
of 8d changes locally, and for the purpose of diagnosis, it 
is desirable that the phase change on the order of 10~ 6 radian 
of 10 urn (= 10" 5 m) is detected. Therefore, the right side 
of Eq. (12) is expressed as follows. 
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d 2 d6 ^ 10~ 6 d<3 
dx 2 (10~ 5 ) 2 



(13) 



Furthermore, substituting -(l/l) x (dx/dz) > 10~ 2 into 
the left side of Eq. (12) and arranging, it is expressed as 
follows . 

1(T 2 ^ d<5 (14) 

5 10~ 10 

d<5 ^ 10" 6 • • • (15) 

Referring to Fig. 4 again, d8 > 10~ 6 is held in the case 
where the X-ray energy is equal or less than about 30 keV. 
Thus, it is desired that the imaging is performed by using 

10 the X-ray having energy from 16 keV to 30keV. 

Next, referring to Figs . 1, 5, and 6, a phase information 
restoring method according to the first embodiment of the 
present invention will be described. Fig. 5 is a flowchart 
showing the phase information restoring method according to 

15 the first embodiment of the present invention. In the 
embodiment, a visible image is constructed by using the 
phase-contrast method based on detection data representing 
two diffraction fringe images taken at the same imaging 
distance by using X-rays having different wavelengths X 0 and 

20 A*. Note that the X-ray having the wavelength X 0 or X ± 
represents a highly monochromatic X-ray with the wavelength 
X 0 or Xi as a center wavelength, and may not strictly be an 
X-ray having single wavelength of X 0 or X x . However , the center 
value of the energy of the X-ray having the wavelength X 0 or 

25 A* is assumed to satisfy the condition from 16 keV to 30 keV. 
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First, at step S10 , X-ray imaging operation is performed. 
That is, as shown in Fig. 6, X-ray imaging is performed by 
positioning the sensor 103 in the position where the imaging 
distance is R and applying an X-ray having wavelength X 0 to 
the object 10. Then, X-ray imaging is performed by applying 
an X-ray having wavelength Xi to the object 10 . Thereby, the 
detection data representing diffraction fringe images are 
obtained . 

By the X-ray imaging at step S10, detection data I (r x , 
R, X 0 ) and I(r x , R, X x ) are sequentially inputted to the image 
constructing unit 200 . Here, the detection data I (r x , R, X Q ) 
represents intensity of the diffraction X-ray having the 
wavelength X 0 in the position r x = (x, y) on a plane at the 
imaging distance of R. Similarly, the detection data I(r x , 
R, X x ) represents intensity of diffraction X-ray having the 
wavelength A* in the position r x = (x, y) on the plane at the 
imaging distance of R . These detection data are sequentially 
stored in the storage unit 201 of the image constructing unit 
200. 

Next, at steps Sll and S12 , the image constructing unit 

200 restores phase <|»(r x , 0, X 0 ) of the X-ray just after 

transmitted through the object on the basis of the detection 

data I(r x , R, X 0 ) and I(r x , R, Xi) stored in the storage unit 
201. 

First, at step Sll, the Laplacian processing unit 202 
obtains a Laplacian f <r x , 0, A, 0 ) = V 2 <|>(r x , 0, Xo) of the phase 
<|>(r x , 0, X 0 ) by using the following equation (16) . 
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f(r ± ,0,A 0 ) = V 2 0(r ± ,O,A o ) = - -^L (ag 0 -a" 2 gi ) .,- (16) 
Where, 

go = ln{I(r_L, R, Xq) } ... (17) 
gi = ln{l(ru R, Ai) } ... (18) 

AX = Xi-X.0, <t = A-iAo 

Accordingly, g 0 and g x are obtained by substituting the 
detection data I(r x , R, Ao) and I{x Xf R, Ai) into Eqs . (17) 
and (18) , respectively, and further, by substituting g 0 and 
gi into Eq. (16) , the Laplacian f (rj., 0, A 0 ) of the phase can 
be obtained. 

Furthermore, at step S12, the inverse Laplacian 
processing unit 203 obtains the phase <j>(r x , 0, A 0 ) by performing 
inverse Laplacian computation on the Laplacian f (r x , 0, Jl 0 ) 
= V 2 <J>(r_L, 0, A, 0 ) of the phase obtained at step Sll . 

Here, the inverse Laplacian computation will be 
described in detail. A Fourier transform on f (r ± , 0, Xq) is 
expressed by the following equation (19) . 
F[f (r ± ,0, A 0 )] = F[V 2 0 (r ± ,0, A 0 )] 

= -4tt 2 (u 2 +v 2 ) [<*> (r ± ,0, A 0 )] •••(19) 
Where, F[ ] represents a Fourier transform and u and v are 
spatial frequencies that correspond to x and y. 

Hereby, the phase <|)(r x , 0, A 0 ) is expressed as Eq. (20) . 

*(rx.0,A o) = P-1[. „„.,,.,,,„] ... <ao , 

Where, F x [ ] represents an inverse Fourier transform. 
Using Eq. (20) , the inverse Laplacian computation can be 
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performed. That is, the restored phase <|>(r]_, 0, Jt 0 ) can be 
obtained by performing the Fourier transform on f (r ± , 0, Jt 0 ) , 
multiplying by { -4n 2 (u 2 + v 2 ) } - 1 and then performing the inverse 
Fourier transform thereon. 

Here, a value of {-4tc 2 (u 2 + v 2 ) J" 1 may be calculated in 
advance within the range where |u| and |v| are not larger 
than a predetermined value so that the calculated value can 
be used when the computation expressedbyEq. (20) isperformed. 
That is, in the case where the predetermined value "const" 
is set, for |u|, |v| < const, the value of the following 
expression is used in Eq. (20) . 

{-4tc 2 (u 2 + v 2 )}"^ (the value calculated in advance) 
On the other hand, for |u| , |v| > const, the value of the 
following expression is used in Eq. (20) . 

-47t 2 (u 2 + v 2 ) } _1 = 0 
Thereby, the inverse Laplacian computation can be performed 
at high speed. 

Next, atstepS13, the image processing unit 204 generates 
image data on the basis of the restored phase 4>(r x , 0, X 0 ) . 
That is, the image processing unit 204 converts the phase 
<|>(rx, 0, X 0 ) in each pixel into data representing brightness , 
and performs necessary image processing such as tone 
processing and interpolation processing. 

Then, as far as necessary, at step S14, the display unit 
300 displays a visible image on the basis of the image data 
on a display, and at step S15, the output unit 400 prints 
it out on a film, etc. 
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E. 

10 



In the embodiment, an X-ray is used when performing 
imaging on an object. However, not limited to the X-ray but 
any radiation that can form diffraction images by being 
transmitted through the object and satisfies the energy 
condition from 16 keV to 30 keV can be used. For example, 
corpuscular beams , etc . including an electron ray can be cited . 
Further, in the embodiment, phase is restored by using two 
X-rays having different energy, however, as described in T. 

Gureyev et ml. "Quantitative In-Line Phase-Contrast 
Imaging with Multienergy X Rays", the phase may be restored 
by using three X-rays having different energy. 

Furthermore, in the embodiment, a synchrotron radiation 
source is used when imaging is performed on an object, a 
radiation source generating beams other than synchrotron 
L5 radiation may be used. For example, the 

electron-storage-type high-brightness hard X-ray generator, 
which has been developed in the Ritsumeikan University, can 
generate X-ray having as high brightness and directionality 
as synchrotron radiation despite of its desktop size . X-ray 
generated by this generator has coherence, and even though 
the X-ray has plural wavelengths , they can be monochromatized 
by combining with monochromatizing crystal. Moreover, the 
ray source developed by The Femtosecond Technology Research 
Association (FESTA) generates ultra short pulse 
high-brightness X-ray on the basis of a principle of backward 
Compton scatter. This ray source is compact and portable, 
and can generate X-ray having not only coherence but also 
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high directionality and monochromaticity . Note that, in the 
case where a point radiation source is used as the radiation 
source , it is desirable to perform correction on the detection 
data obtained by X-ray imaging in view of an enlargement ratio 
before performing data processing in the image constructing 
unit. 

Next, a modified example of the phase information 
restoring apparatus according to the first embodiment of the 
present invention will be described by referring to Fig. 7. 
The phase information restoring apparatus shown in Fig. 7 
has an imaging unit 110 and a reading unit 500. Other 
constitution is the same as that of the phase information 
restoring apparatus shown in Fig. 1. 

In the imaging unit 110, as a screen used for recording 
image information, a photo stimulable phosphor sheet 
(recording sheet) is used instead of the sensor 103 in the 
imaging unit 100 shown in Fig. 2. 

The photo stimulable phosphor (storage phosphor) is a 
material that, when applied with radiation, a part of the 
radiation energy is stored therein, and when applied with 
excitation light such as visible light afterward, stimulated 
fluorescence is generated depending on the stored energy. 
When a radiation image of an object such as a human body, 
etc . is taken and recorded on the sheet that is applied with 
the photo stimulable phosphor, and scanned by the excitation 
light such as laser light, stimulated fluorescent light is 
generated. Therefore, detection data can be obtained by 
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reading the light photoelectrically . After the detection 
data is appropriately processed, the radiation image can be 
displayed as a visible image by being outputted to a display 
such as a CRT or printed out on a film by a laser printer, 
etc. 

The reading unit 500 shown in Fig. 7 is used for reading 
the radiation image recorded on the recording sheet. Here, 
referring to Fig . 8 , constitution and operation of the reading 
unit 500 will be described. The recording sheet 50 on which 
image information has been recorded is set in a predetermined 
position of the reading unit 500. The recording sheet 50 is 
carried in a direction of a Y axis by a sheet carrying means 
52 driven by a motor 51 . On the other hand, a beam LI emitted 
from the laser source 53 is reflected and deflected by a 
rotating multifaceted mirror 55 that is driven by a motor 
54 and rotating at high speed in a direction of an arrow, 
and passes through a convergent lens 56. Then, the beam LI 
has its optical path changed by the mirror 57 and scans the 
recording sheet 50 in a direction of an X axis . By this scanning, 
excitation light L2 is applied to the recording sheet 50, 
and stimulated fluorescent light L3 having quantity depending 
to the stored and recorded radiation image information is 
emitted from the applied part. The stimulated fluorescent 
light L3 is guided by an optical guide 58 and photoelectrically 
detected by a photomultiplier 59. An analogue signal 
outputted from the photomultiplier 59 is amplified by an 
amplif ier 60 and digitized by an A/D converter 61 . The 
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detection data outputted from the A/D converter 61 is inputted 
to the image constructing unit 200. 

Image information representing plural interference 
fringe images obtained in different energy of X-rays can be 
obtained by performing radiation imaging while changing the 
energy of the emitted X-rays and using plural recording sheets 
in the imaging unit 110, and reading image information from 
the respective recording sheets in the reading unit 500 . The 
image constructing unit 200 performs phase restoration on 
the basis of the detection data and generates image data. 
The processing in the image constructing unit 200 is the same 
as that described by referring to Fig. 5. 

Next, a phase information restoring apparatus according 
to a second embodiment of the present invention will be 
described. Fig. 9 is a block diagram showing a constitution 
of the phase information restoring apparatus according to 
the second embodiment of the present invention. 

As shown in Fig. 9, this phase information restoring 
apparatus has an imaging unit 120 for outputting detection 
data that represents radiation image information on an object 
by applying X-rays to the object, an image constructing unit 
210 for generating image data by restoring phase information 
on the basis of the detection data. Other constitution is 
the same as that in Fig. 1. 

Fig . 10 is a diagram showing a constitution of the imaging 
unit 120. The imaging unit 120 includes a radiation source 
101 and a sensor 103 . The X-ray generated from the radiation 
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source 101 is transmitted through the object 10 and enters 
a sensor 103 to generate diffraction fringes. 

The sensor 103 is held by a holding portion 116. The 
holding portion 116 is movably supported on a rail 117. The 
position of the holding portion 116 is controlled by a control 
unit, which will be described later, of the image constructing 
unit 210, and a distance between the object 10 and the sensor 
103 is changed under the control of the control unit. 

Further, the imaging unit 120 includes an amplifier 104 
and an A/D converter 105. The amplifier 104 amplifies 
detection signal outputted from the sensor 103. The A/D 
converter 105 converts the detection signal amplified by the 
amplifier 104 into a digital signal (referred to as "image 
signal" or "detection data") , and outputs the detection data 
to the image constructing unit 210. 

Referring to Fig. 9 again, the image constructing unit 
210 has a storage unit 211 for temporarily storing the detection 
data outputted from the imaging unit 120, a differential 
processing unit 212 for obtaining a differential coefficient 
between detection data at different imaging distances, a 
Laplacian processing unit 213 for calculating a value that 
corresponds to a Laplacian of phase, an inverse Laplacian 
processing unit 214 for performing inverse Laplacian 
computation for performing phase restoration, an image 
processing unit 215 for generating image data on the basis 
of the phase information in a position of the sensor outputted 
from the inverse Laplacian processing unit 214 , and a control 
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unit 216 for controlling the respective units 211-215 and 
the imaging distance in the imaging unit 120. The image 
constructing unit 210 maybe configured with a digital circuit 
or software and a CPU. With a CPU, the control unit 216 
including the CPU processes the detection data on the basis 
of a phase information restoring program recorded on a 
recording medium 217. As the recording medium 2 1 7 , a flexible 
disk, a hard disk, an MO, an MT, a RAM, a CD-ROM, a DVD-ROM, 
etc. are applicable. 

Next, a phase information restoring method according 
to the second embodiment of the present invention will be 
described by referring to Figs . 9-11 . Fig. 11 is a flowchart 
showing the phase information restoring method according to 
the second embodiment of the present invention. In the 
embodiment, a visible image is produced by using image 
information that represents two diffraction fringe images 
taken at different imaging distances z 0 and x x by using X-rays 
having the same wavelength X . Note that the X-ray having the 
wavelength X indicates a highly monochromatic X-ray with the 
wavelength A. as a center wavelength, and may not strictly be 
an X-ray having single wavelength of X. Note that the center 
value of the energy of the X-ray having the wavelength X 
satisfies the condition from 16 keV to 30 keV. 

First, at step S20, the wavelength of the X-ray of the 
X-ray source is set to X, and X-ray imaging is performed while 
changing the position of the sensor 103. That is, the X-ray 
imaging is performed by first positioning the sensor 103 in 
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the position where the imaging distance is z 0 and applying 
an X-ray to the object 10 as shown in Fig. 10 . Then, the X-ray 
imaging is similarly performed after moving the sensor 103 
to the position where the imaging distance is z lm Thereby, 
the image information representing diffraction fringe images 
are obtained. 

By the X-ray imaging at step S20, detection data I(r x , 
z 0 , X) and I(r ± , z lf X) are sequentially inputted to the image 
constructing unit 210. Here, the detection data I(r ± , z 0/ 
X) represents intensity of the diffraction X-ray in the 
position r x = (x, y) on a plane at the imaging distance of 
z 0 . Similarly, the detection data I(r ± , z lf X) represents 
intensity of diffraction X-ray in the position r j_ = (x, y) 
on a plane at the imaging distance of zi. These detection 
data are sequentially stored in the storage unit 211 of the 
image constructing unit 210. 

Next, at steps S21 to S23, the image constructing unit 
210 restores phase in the position of the sensor on the basis 
of the detection data stored in the storage unit 211. 

First, at step S21, the differential processing unit 
212 obtains a differential coefficient between detection data 
I (r x , zi, X) anddetectiondatal (r ± , z 0 , X) by using the following 
equation (21) . 

dl(r ± ,z Q , A) l( r± , Zl , A)-I(r ± ,z 0 , A) 

Tz = — •" (21) 

Then, at step S22, the Laplacian processing unit 213 
obtains a Laplacian f (r x , z, X) = V 2 <|>(r x , z, X) of the phase 
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on the basis of the differential coefficient obtained at step 
S21 and the detection data stored in the storage unit 201 
by using the following equation (22) . 

f(r z A ) - K 3Mr ± ,z 0 , A ) 

1 l r ± ' Z * A ' — ~ " * * ' (22) 

I(r ± ,z 0 , A) dz 
In Eq. (22) , the differential coefficient is divided by the 
detection data I(r_i_, z 0 , X) at the shorter imaging distance. 
However, it may be divided by the detection data I(r X/ z x , 
X) at the longer imaging distance . Further , it may be divided 
by the detection data subjected to LPF (low pass filter) 
processing. 

Furthermore, at step S23, the inverse Laplacian 
processing unit 214 obtains the phase § (r ± , z, X) byperforming 
inverse Laplacian computation on the Laplacian f (r x , z, X) 
= V 2 <t>(r_i_, z, X) of the phase obtained at step S22 . That is, 
the following equation (23) that represents the phase <|>(r ± , 
z, X) is obtained by performing Fourier transform on the 
Laplacian f (r x , z, X) of the phase and transforming it. 

<Mr x ,z,A> = F-[-_L_ F[f(r± , 2 , A)] ] ... (23) 

By using Eq. (23) , the inverse Laplacian computation can be 
performed. By the way, the detailed method of inverse 
Laplacian computation is the same as that described in the 
first embodiment of the present invention. 

Next, atstepS24, the image processing unit 215 generates 
image data on the basis of the restored phase 4>(rj_, z, X) . 
That is, the image processing unit 215 converts the phase 
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<|>(rx, z, A.) in each pixel into data representing brightness, 

and performs necessary image processing such as tone 

processing and interpolation processing. 

Then, as far as necessary, at step S25, the display unit 

300 displays a visible image on the basis of the image data 

on a display, and at step S26, the output unit 400 prints 

it out on a film, etc. 

Next, a phase information restoring apparatus according 

to a third embodiment of the present invention will be described . 
Fig. 12 is a block diagram showing a constitution of the phase 
information restoring apparatus according to the third 
embodiment of the present invention . This phase information 
restoring apparatus has an imaging unit 130 and an image 
constructing unit 220 . Other constitution is the same as that 
of the phase information restoring apparatus shown in Fig. 
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Fig . 13 is a diagram showing a constitution of the imaging 
unit 130. The imaging unit 130 includes a sensor 131. The 
sensor 131 is used as a screen for detecting the incident 
X-ray which produces diffraction fringes, and outputs 
detection signals representing intensity of diffracted X-ray 
entering the respective positions of the sensor 131 . As the 
sensor 131, a two-dimensional sensor such as a CCD (charge 
coupled device) is used, which has a plurality of detecting 
elements that convert intensity of the incident X-ray into 
electric signals and output the signals. 

Further, the imaging unit 130 also includes a sensor 
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driving unit 132 . The sensor driving unit 132 changes the 
distance between the object 10 and the sensor 131 by driving 
a holding portion 116 under control of a control unit 226 
of the image constructing unit 220, which will be described 
later . Other constitution of the imaging unit 130 is the same 
as that of the imaging unit shown in Fig. 10. 

Referring to Fig. 12 again, the image constructing unit 
220 has a storage unit 221 for temporarily storing the detection 
data outputted from the imaging unit 130, a differential 
processing unit 222 for obtaining a differential coefficient 
between detection data at different imaging distances, a 
Laplacian processing unit 223 for calculating a value that 
corresponds to a Laplacian of phase, an inverse Laplacian 
processing unit 224 for performing inverse Laplacian 
computation for performing phase restoration, an image 
processing unit 225 for generating image data on the basis 
of the restored phase information, and a control unit 226 
for controlling the respective units 221-225 and the imaging 
distance in the imaging unit 130 . The image constructing unit 
220 may be configured with a digital circuit or software and 
a CPU. In the latter case, the control unit 226 including 
the CPU processes the detection data on the basis of a phase 
information restoring program recorded on a recording medium 
227. As the recording medium 227, a flexible disk, a hard 
disk, an MO, an MT, a RAM, a CD-ROM, a DVD-ROM, etc. are 
applicable . 

Next, referring to Fig. 14, detection signals detected 
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in a phase information restoring method according to the 
embodiment will be described. As shown in Fig. 14, in the 
embodiment, a combination of a fluorescent screen 133 and 
a CCD camera 134 is used as the sensor 131. That is, the 
fluorescence that is generated by applying the X-ray to the 
fluorescent screen 133 is taken by using the CCD camera 134 
that is disposed oppositely to the X-ray source 101 relative 
to the fluorescent screen 133, thereby an X-ray is detected. 
In the embodiment, the X-ray source 101 having a finite focal 
spot size is used, and the CCD camera 134 with pixels in which 
longitudinal and lateral sizes of the pixels are A x and A Y , 
respectively . 

Here, a blur is produced in the image obtained by the 
X-ray generated from the X-ray source 101 having the finite 
focal spot size a. Assuming that the distance between the 
object 10 and the X-ray source 101 is R, the imaging distance 
between the object 10 and the sensor 131 is z, and the focal 
spot sizes of the X-ray source in directions of an x axis 
and an y axis are cr x and ct y , respectively, blur function f (u, 
v) representing the blur of the image assumes normal 
distribution expressed by the following equation. 

f(u,v) = exp[-ya'u 2 -Y^v 2 ] •••(24) 

Where, u and v are spatial frequency components in the 
directions of the x axis and the y axis in the sensor 131, 
respectively, and a x = 2toj x z/R, a Y = 27iCT y z/R. 

Fig . 15 shows the blur function as a function of spatial 
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frequency . As shown in Fig . 15 , the value of the blur function 
in the high spatial f requency region takes approximately zero . 
Therefore, even if the pixel size of the detecting elements 
that constitute the sensor is made smaller in order to obtain 
high-definition images, the detection signals of the images 
become few under the influence of blurs in the high spatial 
frequency region, thereby only noise is detected. 

On this account, in order to reduce the influence of 
noise , the highest spatial frequency that the sensor can detect 
may be made lower than the spatial frequency at which the 
detection signals become few due to blurs . That is , assuming 
that the standard deviation in the blur function is sigma 
( <T) , the highest spatial frequency that the sensor can detect 
maybe within the range of 3 sigma including 99% of the detection 
signals . 

FromEq. (24) , the triple sigma (3(7) range is represented 
by u < 3/a x , v < 3/a Y , and the highest spatial frequency relative 
to the pixel sizes A x , A y of the sensor is representedby 1/ (2A X ) , 
1/ (2A Y ) according to the sampling theorem. Therefore, it is 
understood that the pixel sizes A x , A Y may satisfy the following 
equation in order to reduce the influence of noise. 
1 < 3 1 ^ 3 

2A X ^' "^7-17 " (25) 

By arranging the Eq. (25) , the following equation (26) 
that shows the condition of the pixel sizes required for 
reducing the influence of noise can be obtained. 

A x > a x /6, A Y > a Y /6 . . . (26) 
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For example, in the case where z = 1.0m, R = 5.0 m, a x = a Y 
= 40 nm, the pixel sizes may be equal or more than 8.4 pm. 

Next, referring to Figs. 12, 13, and 16, the phase 
information restoring method according to the third 
embodiment of the present invention will be described. Fig. 
16 is a flowchart showing the phase information restoring 
method according to the third embodiment of the present 
invention. In the embodiment, diffraction images are taken 
twice while changing the imaging distance , and a visible image 
is constructed on the basis of the detection data representing 
image information on these diffraction images by using the 
phase-contrast method. 

First , at step S30 , X-ray imaging operation is performed . 
That is, as shown in Fig. 13, the object is positioned in 
a position where the distance from the X-ray source 101 is 
R, and the sensor driving unit 132 positions the sensor 131 
in a position where the imaging distance is z x under the control 
of the control unit 226, and an X-ray is applied to the object 
10 , thereby X-ray imaging is performed . Similarly, the X-ray 
imaging is performedby positioning the sensor 131 inaposition 
where the imaging distance is z 2 . 

Here, assuming that the focal spot sizes of the X-ray 
source in directions of X, Y are ct x , a Y , both the longitudinal 
and lateral pixel sizes of the detecting element of the sensor 
131 are A, the larger one of the focal spot sizes ct x , a y is 
a, and the longer one of the imaging distances is z , the sensor 
driving unit 132 positions the sensor 131 using the detecting 
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element having the pixel size that satisfies A > 7taz/3R so 
that Eq. (26) may be satisfied. 

By the X-ray imaging at step S30, the detection data 
I(x, y, zi) and I(x, y, z 2 ) that represent intensity of 
diffraction X-ray entering the pixel (x, y) on planes at the 
imaging distance of z x and z 2 , respectively, are sequentially 
inputted to the image constructing unit 220, and stored in 
the storage unit 221. These detection data represent 
diffraction fringe image information on the planes at the 
respective imaging distances . 

Next, at steps S31-S33, the image constructing unit 220 
restores phase <|>(x, y) in the position of the sensor on the 
basis of the detection data I(x, y, Zl ) and I(x, y, z 2 ) stored 
in the storage unit 221. 

First, at step S31, the differential processing unit 
222 obtains a differential coefficient between the detection 
data I(x, y, Zl ) and the detection data I(x, y, z 2 ) using the 
following equation (27) . 

dl(x,y,z) l(x,y,z 2 )-l(x,y, 2l ) 

- = ■ • • (27) 

dz z 2 _ z i 

Then, at step S32, the Laplacian processing unit 223 
obtains Laplacian f (x, y, z) = V 2 (|>(x, y, z) of phase on the 
basis of the differential coefficient obtained at step S31 
and the detection data stored in the storage unit 221 using 
the following equation (28) . 

f (v „ K dl(x, y,z) 

r(x,y,z)- _ •••(28) 

I(x,y, Zl ) dz K a) 
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Here, in Eq. (28) , although the differential coefficient is 
divided by the detection data I (x, y, Zl ) at a shorter imaging 
distance, it may be divided by the detection data I(x, y, 
z 2 ) at a longer imaging distance. 

Further, at step S33, the inverse Laplacian processing 
unit 224 obtains phase <J>(x, y, z) by performing inverse 
Laplacian computation on the Laplacian f (x, y, z) = V 2 <|>(x, 
y , z) of the phase obtained at step S32 . Note that the method 
of the inverse Laplacian computation is the same as that 
described in the first embodiment of the present invention. 

Next, atstepS34, the image processing unit 225 generates 
image data on the basis of the phase <|>(x, y, z) . That is, 
the image processing unit 225 converts the phase <|>(x, y, z) 
in each pixel into data representing brightness and performs 
necessary image processing such as tone processing and 
interpolation processing. 

Then, as far as necessary, at step S35, the display unit 
300 displays a visible image on the basis of the image data 
on a display, and at step S36, the output unit 400 prints 
it out on a film, etc. 

In the embodiment, an X-ray is used when performing 
imaging on an object. However, not limited to the X-ray but 
any radiation that can form diffraction images by being 
transmitted through the object, for example, corpuscular 
beams, etc. including an electron ray can be used. 
Alternatively, similarly to that described in the first 
embodiment of the present invention, a radiation source 
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generating beams other than synchrotron radiation maybe used. 

In addition, in the embodiment, the phase is restored 
by using two pieces of detection data at different imaging 
distances . However, phase may be restored by using equal or 
more than three pieces of detection data at different imaging 
distances . 

Further, when performing radiation imaging, the X-ray 
transmitted through the object lOmaybe enlarged and detected. 
For example, the fluorescence generated from the fluorescent 
screen 133 shown in Fig. 14 is enlarged by an optical system 
so as to be taken by the CCD camera 134. In this case, the 
effective pixel size (pixel size before enlarged) is used 
as A x and Ay instead of the actual pixel size of the CCD. 
Alternatively, as shown in Fig. 17, by using a point source, 
etc. as the X-ray source, the X-ray may be detected in an 
enlarged scale, and also in this case, the effective pixel 
size is similarly used as A x and Ay. 

Next, a phase information restoring apparatus according 
to a fourth embodiment of the present invention will be 
described. Fig. 18 is a block diagram showing a constitution 
of the phase information restoring apparatus according to 
the fourth embodiment of the present invention. 

As shown in Fig. 18, this phase information restoring 
apparatus includes an image constructing unit 230 instead 
of the image constructing unit 220 shown in Fig. 12 , and other 
constitution is the same as that shown in Fig. 12 . The image 
constructing unit 230 is constituted by further adding an 
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LPF (low pass filter) 231 for performing low pass filter 
processing on the detection data outputted from the imaging 
unit 130 under the control of the control unit 226 and 
outputting the data to the storage unit 221 to the image 
constructing unit 220 shown in Fig. 12. In the embodiment, 
the low pass filter processing is performed by digital 
computation in order to suppress or eliminate high spatial 
frequency components . 

Next, a phase information restoring method according 
to the fourth embodiment of the present invention will be 
described by referring to Figs. 13, 18, and 19. Fig. 19 is 
a flowchart showing the phase information restoring method 
according to the fourth embodiment of the present invention . 
In the embodiment, diffraction images are taken twice while 
changing the imaging distance and a visible image is 
constructed on the basis of the detection data representing 
image information on these diffraction images by using the 
phase-contrast method. 

First, atstepS40, X-ray imaging operation is performed. 
That is, as shown in Fig. 13, the object is positioned in 
a position where the distance from the X-ray source 101 is 
R, and the sensor driving unit 132 positions the sensor 131 
in a position where the imaging distance is z ± under the control 
of the control unit 226, an X-ray is applied to the object 
10 , thereby X-ray imaging is performed . Similarly , the X-ray 
imaging is per formed by positioning the sensor 131 inaposition 
where the imaging distance is z 2 . 
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By the X-ray imaging operation in step S40 , the detection 
data I(x, y, and I(x, y, z 2 ) representing intensity of 

the diffraction X-ray entering the pixel (x, y) on the planes 
at respective imaging distances Zi and z 2 are sequentially 
inputted to the image constructing unit 220 . These detection 
data represent diffraction fringe image information on the 
planes at the respective imaging distances. 

Next, at step S41, the LPF 231 performs low pass filter 
processing on the sequentially inputted detection data so 
as to obtain detection data I LPF (x, y, z x ) and I LPP (x, y, z 2 ) 
in which high spatial frequency components of the detection 
data are suppressed or eliminated . The detection data I LPP (x , 
y, zi) andI LPF (x, y, z 2 ) obtained by the LPF 231 are sequentially 
stored in the storage unit 221. 

Here, the region in which the high spatial frequency 
components are suppressed or eliminated by the LPF 231 is 
a region in which there is no signal because of the influence 
of blurs , as shown in Fig . 15 . Therefore , in order to suppress 
or eliminate the spatial frequency components outside the 
triple sigma range, the spatial frequency components that 
satisfy u > 3R/27kt x z, v > 3R/27i;aYZ may be suppressed or 
eliminated. Where, the longitudinal and lateral focal spot 
sizes are a x , a Y , the distance between the X-ray source and 
the object is R, the imaging distance is z, and the spatial 
frequency components of x, y are u, v. 

Next, at steps S42-S44, the image constructing unit 230 
restores phase <|>(x, y) in the position of the sensor on the 
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basis of the detection data Ii lPF (x, y, 2l ) and I LPF (x, y, z 2 ) 
stored in the storage unit 221. 

First, at step S42, the differential processing unit 
222 obtains a differential coefficient between the detection 
data I LPF (x, y, Zi) and the detection data I LPF (x, y, z 2 ) using 
the following equation (29) . 

dl(x,y,z) I L p P (x /y/ z 2 ) - I LPF (x,y, 2l ) 

T* — = — ■■'< 29 > 

Then, at step S43, the Laplacian processing unit 223 
obtains Laplacian f (x, y, z) = V 2 <|>(x, y, z) of phase on the 
basis of the differential coefficient obtained at step S42 
and the detection data stored in the storage unit 221 by using 
the following equation (30) . 

f(x,y,z)= ■••(30) 

i lpf( x /Y^ z i) dz 

Here, in Eq. (30) , although the differential coefficient is 

divided by the LPF processed detection data I L pf (x,y , zi) at 

the shorter imaging distance, it may be divided by the LPF 

processed detection data I LPF ( X/ y , z 2 ) at the longer imaging 

distance . 

Further, at step S44, the inverse Laplacian processing 
unit 224 obtains phase <|>(x, y, z) by performing inverse 
Laplacian computation on the Laplacian f (x, y, z) = V 2 <|>(x, 
y , z) of the phase obtained at step S43 . Note that the method 
of the inverse Laplacian computation is the same as that 
described in the first embodiment of the present invention. 

Next, atstepS45, the image processing unit 225 generates 
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image data on the basis of the phase <|>(x, y, z) . That is, 
the image processing unit 225 converts the phase <|>(x, y, z) 
in each pixel into data representing brightness and performs 
necessary image processing such as tone processing and 
interpolation processing. 

Then, as far as necessary, at step S46, the display unit 
300 displays a visible image on the basis of the image data 
on a display, and at step S47, the output unit 400 prints 
it out on a film, etc. 

Next, a modified example of the phase information 
restoring apparatus according to the second to fourth 
embodiments of the present invention will be described. In 
the phase information restoring apparatus according to the 
second to fourth embodiments of the present invention, in 
stead of the sensor, a photo stimulable phosphor sheet 
(recording sheet) can be used as the screen used for recording 
image information similarly to that in the first embodiment. 

Fig . 20 shows a modified example of the phase information 
restoring apparatus according to the fourth embodiment of 
the present invention. This phase information restoring 
apparatus has an imaging unit 140 and a reading unit 500. 
Other constitution is the same as that of the phase information 
restoring apparatus shown in Fig. 18. 

In the imaging unit 140, instead of the sensor 131 shown 
in Fig. 13, a photo stimulable phosphor sheet (recording sheet) 
is used. Further, the constitution and the operation of the 
reading unit 500 are similar to those described by referring 
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to Fig . 8 . 

In the imaging unit 140, radiation imaging is performed 
by changing the imaging distance and using plural recording 
sheets. Further, in the reading unit 500, image information 
is read from the respective recording sheets, thereby 
detection data representing plural diffraction fringe images 
obtained at different imaging distances can be obtained . The 
image constructing unit 230 restores phase on the basis of 
this detection data and generates image data . The processing 
in the image constructing unit 230 is the same as that described 
by using Fig. 19. 

As described above, according to the present invention, 
when constructing the radiation image of a living organism 
such as a human body by the phase -contrast method, the 
estimation accuracy of phase can be improvedby using radiation 
with energy of high transmittance from 16 keV to 30 keV. In 
addition, according to the present invention, when 
constructing the radiation image of a living organism such 
as a human body by the phase-contrast method, the influence 
of noise can be reduced without increasing irradiation amount 
of the X-ray by eliminating the detection signals in the region 
where there is no signal because of blurs. 
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